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ABSTRACT 

It has been proposed that very massive stars (VMSs) dominated heavy ele- 
ment production until a "metallicity" threshold corresponding to [Fe/H] ~ —3 
was reached. This results in a prompt (P) inventory of elements, the abun- 
dances of which were determined from observations of Galactic halo stars with 
[Fe/H] w —3. We calculate Qp(E) from the P inventory for a large number 
of elements in the IGM. Using the available data on f2(E; on ) for C IV, VI, 
and Si IV in Lya systems, we find that the ionization fractions calculated from 
f2(E ion )/f2p(E) are, within reasonable uncertainties, compatible with the values 
estimated from ionization models for Lya systems. This agreement appears to 
hold from z ~ 0.09 to ~ 4.6, indicating that the bulk of the baryonic mat- 
ter remains dispersed with a fixed chemical composition. We conclude that the 
P inventory was established in the epoch prior to z ~ 4.6. The dispersal of 
processed baryonic matter to the general IGM is considered to be the result of 
energetic VMS explosions that disrupted baryonic aggregates until the "metal- 
licity" threshold was reached to permit normal astration. The formation of most 
galaxies is considered to have occurred subsequent to the achievement of this 
metallicity threshold in the IGM. 

Subject headings: early universe — intergalactic medium — quasars: absorption 
lines 



1. Introduction 

We present predictions for the abundances of a large number of elements in the inter- 
galactic medium (IGM), especially the Lya systems, based on considerations of a prompt 
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(P) inventory that is dominated by the contributions from very massive (> 100 M & ) stars 
( VMSs) . A comparison will be made between the predictions and the available observational 
data. The model of the P inventory (Wasserburg & Qian 2000a; Qian & Wasserburg 2001a; 
Qian & Wasserburg 2002) claims that, prior to the achievement of the condition correspond- 
ing to a metallicity of [Fe/H] ps —3 in the IGM, the predominant mechanism of astration is 
by formation of VMSs. This condition is identified based on the sharp rise in the abundances 
of the heavy r-process elements such as Ba and Eu in Galactic halo stars with [Fe/H] ps —3 
(e.g., McWilliam et al. 1995; Burris et al. 2000). This sharp rise is considered to represent 
the occurrence of Type II supernovae (SNe II) that result from evolution of normal stars 
with masses of ~ (10-60) M Q in the absence of VMSs. Further support for a threshold at 
[Fe/H] ps —3 for normal astration comes from the observation that all damped Lja systems 
with redshift z ~ 1.5-4.5 have a minimum [Fe/H] of ps —2.7 (e.g., Prochaska & Wolfe 2000; 
Prochaska, Gawiser, & Wolfe 2001). The metallicity of [Fe/H] ps —3 is considered to rep- 
resent the critical condition for sufficient cooling and fragmentation of gas clouds to permit 
major formation of normal stars with masses of ~ (1-60) M (Wasserburg & Qian 2000a). 
This is supported by ab initio calculations of Bromm et al. (2001), which showed that a 
critical metallicity of ~ 5 x 10~ 4 the solar value (corresponding to [Fe/H] ~ —3) is required 
for significant formation of low-mass protostellar aggregates. 

It was found that the observed elemental abundances in stars with —3 < [Fe/H] < — 1 
could be quantitatively explained by the contributions from: (1) the P inventory at [Fe/H] ps 
—3 for all the elements up to Zr; (2) a high-frequency type of SNe II [SNe 11(H)] that are 
responsible for the heavy r-process elements and some light r-process elements below Ba but 
produce very little of the elements from N to Zn including Fe; and (3) a low-frequency type 
of SNe II [SNe II(£)] that are responsible for the bulk of the light r-process elements and 
the part of the elements from N to Zn attributable to SNe II (Qian & Wasserburg 2001b, 
2002). It was further found by reducing the observational data on metal-poor Galactic halo 
stars that the yield patterns of SNe II (L) and VMS are almost the same for all the elements 
below Sr except for some small but significant shifts at Na, Al, V, Cr, Mn, and Co. Thus, 
the abundance ratios of the so-called a elements relative to Fe such as Mg/Fe, Si/Fe, Ca/Fe, 
and Ti/Fe remain essentially constant at —4 < [Fe/H] < —1 although the dominant source 
of the a elements and Fe changes from VMSs in the regime of [Fe/H] < —3 to SNe ll(L) in 
the regime of —3 < [Fe/H] < — 1. These ratios are all greater than the corresponding solar 
values by a factor of ~ 3 because ~ 2/3 of the solar Fe inventory is from Type la Supernovae, 
which later added the Fe group elements but little of the a elements at [Fe/H] > — 1. This 
simple explanation should serve to clarify the so-called "overproduction" of what are referred 
to as the a elements at low metallicities. 

It was considered that the P inventory at [Fe/H] ps —3 is dominated by the contributions 
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from VMSs with small contributions from SNe II. The rare occurrence of Galactic halo stars 
with —4 < [Fe/H] < —3, which must have masses of ~ 1 M Q in order to survive to the present 
epoch, is considered to represent concomitant formation of some "normal" stars in the regime 
where VMSs dominated astration (Qian & Wasserburg 2002). It was also suggested that 
VMSs were the source for reionization of the IGM (e.g., Qian & Wasserburg 2001b). A recent 
study by Oh et al. (2001) showed that in producing a "metal" inventory corresponding to 
[Fe/H] ~ —3 (or [Si/H] ~ —2.5), VMSs also provide ~ 1-10 ionizing photons for H and He. 
This was estimated from the VMS model yields calculated by Heger & Woosley (2002) and 
the VMS emission spectra calculated by Bromm, Kudritzki, & Loeb (2001). It appears that 
by processing ~ 10~ 5 -10~ 4 of the baryonic matter, VMSs are capable of providing the P 
inventory and also of reionizing the IGM prior to the formation of most galaxies. 

Further justification for the role of VMSs is found from the very low Ba abundances at 
[Fe/H] < —3. This requires that any Ba produced by SNe II at [Fe/H] < —3 had to be mixed 
in a dilution mass of ~ (10 6 -10 7 ) M Q , to be compared with the dilution mass of ~ 3 x 10 4 M® 
for SNe II at [Fe/H] > —3. It is considered that the severely disruptive VMS activities at 
[Fe/H] < —3 would not often permit the preservation of baryonic aggregates with masses 
of < (10 6 -10 7 ) M Q and that the P inventory at [Fe/H] « —3 should represent the chemical 
composition of the general IGM (Qian & Wasserburg 2002). Elemental abundances in the 
IGM are best measured from absorption lines produced by regions of enhanced density such 
as Lya systems that are illuminated by Quasi-Stellar Objects (QSOs). For a given QSO, 
absorption lines occur for narrow ranges in redshift (below z of the QSO) and the absorbing 
systems represent concentrations of dispersed baryonic matter between the QSO and the 
observer. Extensive studies have determined the abundances of C IV, O VI, and Si IV in 
Lya systems with a wide range in redshift from z ~ 0.09 up to z ~ 5.3. In the following, we 
will present the abundances of a large number of elements in the IGM based on the model 
of the P inventory (§2) and compare these results with the available data on Lya systems 
(§3). The cosmological implications of our model will be discussed (§4). 

2. Elemental Abundances in Ly« Systems as Determined from the P 



The abundances for the P inventory are given in Table 1 in the standard spectroscopic 
notation as 



where (E/H)p is the number ratio of an element E relative to hydrogen for the P inventory. 
The logep(E) values in Table 1 are directly obtained from observational data on Galactic 



Inventory 




(1) 
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halo stars with [Fe/H] —3 and correspond to those given in Qian & Wasserburg (2002). 
The only element for which there may be a serious uncertainty is O. The logep(O) value 
assumed here is based on observations of HD 115444 and HD 122563 (Westin et al. 2000). 
Other extensive studies (e.g., Israelian, Garcia Lopez, & Rebolo 1998; Boesgaard et al. 1999) 
indicate that logep(O) may be up to ~ 0.3 dex higher than the value given in Table 1. The 
precise value of logep(O) is an important issue that remains to be resolved. 

We consider that essentially all of the baryonic matter was ionized and dispersed in 
the IGM when VMS activities ceased and that the chemical composition of this baryonic 
matter was represented by the P inventory. The mechanism for maintaining the ionized 
state of the IGM will not be discussed here. We assume that subsequent to the acquisition 
of the P inventory, the bulk of the baryonic matter remains dispersed in the form of Lja 
systems. We further assume that these systems have experienced no further evolution in 
elemental abundances and that they represent regions of enhanced density comoving with 
the expansion of the universe. Under these assumptions, the number ratio of an element E 
relative to H in the Lja systems is the same as that for the P inventory and can be written 
as 

/E\ n P (E)/m E = n P (E)/A E 

\Rj p ft(H)/m H ft(H) ' 1 j 

where typ(E) and f2(H) are the fractions of the critical mass density p cri = 3Hq/(8itG) 
contributed by E and H, respectively, at the present epoch, m E and m H are the atomic 
masses of E and H, respectively, and A E = m E /mn is the atomic mass number of E relative 
to hydrogen. The symbols H and G in the expression of p cr j are the Hubble and gravitational 
constants, respectively. The baryonic contribution to p cvl is found to be ~ 0.02/i~ 2 based 
on measurements of deuterium abundances in QSO absorption systems (e.g., Buries, Nollett, 
& Turner 2001) and the power spectrum of the cosmic microwave background (e.g., de 
Bernardis et al. 2002), where h is H in units of 100 km s" 1 Mpc" 1 . This gives f2(H) 
X u il b ~ 0.0152/i~ 2 , where X H ~ 0.76 is the mass fraction of H. Combining equations (1) 
and (2), we obtain 

ftp(E) = Q(E)A E x io lo s^( E )- 12 w 1.52h- 2 A E x io lo s^( E )- 14 . (3) 

Note that fl P (E) does not explicitly depend on cosmological models except for the numerical 
values of H and fife. The values of fip(E) calculated from the P inventory are given in Table 
1 for h = 0.65. At present, data are available only for several elements such as C, O, and Si 
in Lya systems. 
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3. Comparison with Observational Data 

Absorption lines produced by ions such as C IV, O VI, and Si IV in regions of enhanced 
density such as the Lya systems have been measured in the spectra of QSOs that illuminate 
such regions. The column density of the ions of an element E, E ion , in a redshift interval 

z_ < z < z + is 

"+ n(E„ 



where c is the speed of light, n(E ion , z) is the number density of E ion at redshift z, and H(z) 
is the Hubble parameter at z. Under the assumptions that the bulk of the baryonic matter 
resides in Lyo systems and that these systems experience no evolution other than comoving 
with the expansion of the universe, we obtain 

where n (E ion ) = n(E ion , z)/(l + z) 3 is the number density of E ion at the present epoch. By 
defining (Bahcall & Peebles 1969) 

and AX = X(z + ) — X(z_), equation (5) can be rewritten as 

X(E ion ) = Cn °[f ion) AX. (7) 

-no 

If absorption lines due to E ion have been measured for a statistical sample of Lya systems 
at different redshift along the lines of sight for a number of QSOs, the terms X(E ion ) and 
AX in equation (7) are replaced by the sums Yl X(E ion ) and Yl AX over the sampled 
region. It follows that the fraction of p cri contributed by E ion at the present epoch, f2(E ion ) = 
n (E ion )m E /pcri, is 

Won) = ^^^-e = 9.55 x lO-*fc-U E ^^, (8) 

c Pcri 2^ AA AA 

where the last equality is obtained for Y ^(Ei on ) i n units of cm" 2 . If the model of the P 
inventory is correct, the value of fi(Ei on ) obtained from equation (8) is related to that of 
Qp(E) from equation (3) through the ionization fraction /(E ion ) = f2(E ion )/f2p(E). 

The calculation of f2(E ion ) through AX requires a specific cosmological model. This is 
distinct from the more direct calculation of Q P (E). The choice of a particular model makes 
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little difference for z <C 1, as all models give AX m z + — z_ in this limit. For simplicity in 
considering the data on Lja systems with a wide range in z, we adopt a flat cosmological 
model with Q m + = 1, where Q m and are the fractions of p CT \ contributed by matter 
(baryonic and nonbaryonic) and the cosmological constant, respectively, at the present epoch. 
With this model, [H(z)/H Q ] 2 = Q m (l + z) 3 + Q\ and the expression of X(z) from equation 
(6) is 

x(z) = Un m (i + zf + i- n m - 1] , (9) 

which gives AX ^ z + - z„ for z < 1 and AX « |fi m 1/2 [(1 + z + ) 3 / 2 - (1 + z„f/ 2 ] for 
l + ^>[(l-Q m )/Q m ]V3. 

A comparison of the observed values of f2(E ion ) with the values corresponding to the P 
inventory requires that the same cosmological model and parameters be used. The data on 
Lya systems with a wide range in z as reported in the literature were based on somewhat dif- 
ferent cosmological models and parameters. We have recalculated f2(E ion ) from the published 
data so that the recalculated results correspond to the flat cosmological model adopted here 
with Q m = 0.3, Q\ = 0.7, and h = 0.65. The data on C IV and Si IV are taken from Songaila 
(2001) and those on O VI from Buries & Tytler (1996), Tripp, Savage, & Jenkins (2000), 
Savage et al. (2002), and Simcoe, Sargent, & Rauch (2002; no recalculation necessary). The 
recalculated f2(E ion ) values are given in Table 2. Ideally, the fraction /(E ion ) of the element 
E that is ionized to the particular state observed can be calculated for the Lya systems from 
which fi(Eion) is obtained. A comparison can then be made of Q(E) = r2(E ion )//(E ion ) with 
Q P (E) calculated from the P inventory. As values of /(E ion ) are not always well known, 
we have tabulated the ratio fi(E ion )/fi P (E) for C IV, O VI, and Si IV in Table 2. This 
ratio should correspond to /(Ei on ) if the model of the P inventory is valid. The values of 
fi(E ion )/fi P (E) in Table 2 are in reasonable agreement with f(C IV) - 0.5, f(0 VI) ~ 0.2, 
and /(Si IV) ~ 0.1 as estimated from various ionization models for Lya systems. 

For all the data sets used in the above comparison, there may be significant uncertainties 
in estimates of both -^(Eion) and Yl AX due to sampling problems. The mechanisms of 
ionization are also not well understood. Considering the very different and complex data sets 
used and the extreme simplicity of our model for the P inventory, we conclude that there 
is remarkable agreement between the observations and the predictions for the abundances 
in the IGM. We note that the apparently high values of (O/H)/(O/H) and (C/H)/(C/H) 
reported by Buries & Tytler (1996) are off due to the omission of a factor l/A E in their 
equation for ((z) (Tytler 2002, personal communication). We also note that the values of 
[O/H] used by Tripp et al. (2000) and Savage et al. (2002) in estimating the lower bounds 
to fib may not be appropriate in view of the results presented here. We consider that a 
more detailed comparison of the model for the P inventory and the observations must await 
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the availability of more complete data sets with appropriate estimates of errors including 
those in the term AX. It would be most useful to theoretical studies of Lya systems 
if observational reports included explicit statements on the value of iV(E ion ) and on the 
detailed redshift intervals for calculation of ^ AX with corrections for the obscured part of 
the redshift path. We note that compared with the flat cosmological model with Q m = 0.3 
and Q\ = 0.7 as adopted here, a flat model with fl m = 1 (i.e., with the deceleration parameter 
q = 0.5) gives smaller AX values, and hence larger f2(E ion ) values, by a factor of 1.83 at 
high z. It is possible that the model for the P inventory may be of use in evaluating various 
cosmological models if the ionization mechanism for Lya systems is better understood. 



4. Discussion 

The results in Table 2 strongly support our assumptions that subsequent to the ac- 
quisition of the P inventory, the bulk of baryonic matter remains dispersed in the form of 
Lya systems and that these systems have experienced little or no chemical evolution but 
are simply comoving with the expansion of the universe. It is striking to note that these 
assumptions hold from z ~ 2.5 down to ~ 0.09 based on the O VI data and from z ~ 4.6 
down to ~ 2 based on the C IV and Si IV data. We conclude that the bulk of ionized 
baryonic matter in the IGM has not been greatly diminished from the epoch of reionization 
until the present epoch. 

We have argued that VMSs are responsible for the P inventory and reionization of the 
IGM. An important issue is the timescale over which this is achieved. From the range in 
z where the model appears to describe the abundances observed in Lya systems, we infer 
that the IGM had already acquired the P inventory before z ~ 4.6. This is consistent with 
the following considerations of the cosmological epoch during which VMSs could disperse 
processed baryonic matter to the general IGM. The explosion energies of VMSs are ~ 10 52 - 
10 53 erg as estimated from the available models of Heger & Woosley (2002). Based on 
the cold dark matter model of structure formation, dark matter halos corresponding to 
1 a fluctuations and collapsing at z ~ 4 or those corresponding to 2 a fluctuations and 
collapsing at z ~ 10 provide binding energies of ~ 10 52 -10 53 erg for the baryonic matter in 
their potential wells (see Figure 9 in Barkana & Loeb 2001). Thus, VMS explosions could 
readily unbind baryonic matter that was not in stars from most dark matter halos formed at 
z > 4. We conclude that VMSs dominate both the elemental production and the dispersal of 
processed baryonic matter at z > 4. Compared with the nearly constant values of f2(C IV) 
and f2(Si IV) at z ~ 2-4.6, the apparent decrease in these two quantities from z ~ 4.6 to 
~ 5.3 (see Table 2 here or Table 1 in Songaila 2001) may be a hint for cessation of VMS 
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activities at z ~ 4.6 although this decrease may also be due to the incompleteness of the 
sample of Lya systems in the highest redshift interval. Songaila and Cowie (2002) have 
reported a damped Lya system with [Fe/H] ~ —2.7 and z = 5.3. This indicates that VMS 
activities may have ceased as early as z ~ 5.3. In any case, the bulk of baryonic matter 
should have been enriched with the P inventory, ionized, and dispersed in the IGM when 
the epoch of VMS activities ended. The subsequent reaggregation of baryonic matter into 
galaxies appears to be very inefficient as demonstrated by the prevalence of the P inventory 
in Lya systems with z ~ 0.09-4.6. When galaxies do form, it appears that the times of 
formation vary greatly, as argued by Wasserburg & Qian (2000b) based on the large scatter 
in [Fe/H] for damped Lya systems at a given redshift in the range of z « 1.5-4.5. The 
precise epoch of VMS activities and its influence on structure formation at later times are 
important issues that remain to be addressed. 

We would like to dedicate this paper to the memory of Fritz Zwicky. We thank An- 
toinette Songaila, Todd Tripp, and David Tytler for generously answering many questions 
and for guiding us through their data. This work was supported in part by the Department of 
Energy under grants DE-FG02-87ER40328 and DE-FG02-00ER41149, and by NASA under 
grant NAG5-10293, Caltech Division Contribution 8780(1091). 
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Table 1. fip(E) Values as Determined from the P Inventory 



E 


A E 


logep(E) 


n P (E) 


[E/H] 


(1) 


(2) 


(3) 


(4) 


(5) 


C 


12 


5.11 a 


5.56 x 10~ 8 


-3.45 


N 


14 


6.30 


1.00 x 10~ 6 


-1.75 


O 


16 


6.60 


2.29 x 10~ 6 


-2.33 


Na 


23 


3.34 


1.81 x 10~ 9 


-2.99 


Mg 


24 


5.13 


1.16 x 10~ 7 


-2.45 


Al 


27 


3.12 


1.28 x 10~ 9 


-3.35 


Si 


28 


5.02 


1.05 x 10~ 7 


-2.53 


Ca 


40 


3.75 


8.09 x 10~ 9 


-2.61 


Sc 


45 


0.28 


3.08 x 10~ 12 


-2.82 


Ti 


48 


2.43 


4.65 x 10~ 10 


-2.56 


V 


51 


1.15 


2.59 x lO" 11 


-2.85 


Cr 


52 


2.40 


4.70 x 10~ 10 


-3.27 


Mn 


55 


1.90 


1.57 x 10~ 10 


-3.49 


Fe 


56 


4.51 


6.52 x 10~ 8 


-3.00 


Co 


59 


2.24 


3.69 x lO" 10 


-2.68 


Ni 


58 


3.25 


3.71 x 10~ 9 


-3.00 


Cu 


63 


0.53 


7.68 x 10~ 12 


-3.68 


Zn 


64 


1.86 


1.67 x 10~ 10 


-2.74 


Sr 


88 


0.13 


4.27 x 10~ 12 


-2.77 


Ba 


138 


-1.80 


7.87 x 10" 14 


-3.93 



Note. - Column (2) gives the approximate 
atomic mass number for the element; (3) the P in- 
ventory; (4) the fraction of the critical mass density 
contributed by the element at the present epoch; 
and (5) [E/H] = loge(E) - loge (E). 

a The log ep value for C is obtained from data on 
HD 115444 and HD 122563 (Westin et al. 2000). 
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Table 2. Ionization Fractions for Lya Systems 



E 


fi(E) 


E- 

Moil 


{*) 


f^(Eion) 


n(E ion )/n P (E) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


C 


5.56 x 1CT 8 


C IV 


1.87 


2.05 x 10~ 8 


0.37 








2.27 


3.74 x 10~ 8 


0.67 








2.78 


2.62 x 10~ 8 


0.47 








3.21 


2.46 x 10~ 8 


0.44 








3.75 


3.68 x 10~ 8 


0.66 








4.24 


3.93 x 10~ 8 


0.71 








4.66 


1.35 x 10~ 8 


0.24 








5.29 


4.93 x 10~ 9 


0.09 


o 


2.29 x lCT fi 


O VI 


0.09 


4.6 x 10~ 7 


0.20 








fa 0.24 


1 x 10~ 6 


0.44 








0.9 


> 1.3 x 10~ 7 


> 0.06 








» 2.5 


2.96 x 10~ 7 


0.13 


Si 


1.05 x icr 7 


Si IV 


2.24 


9.04 x 10~ 9 


0.09 








2.81 


1.59 x 10~ 8 


0.15 








3.19 


6.19 x 10~ 9 


0.06 








3.70 


1.59 x 10~ 8 


0.15 








4.28 


1.63 x 10~ 8 


0.16 








4.64 


4.66 x 10~ 9 


0.04 








5.29 


9.31 x lO" 10 


0.01 



Note. — Data on C IV and Si IV from Songaila 2001 and data on O VI 
from Savage et al. 2002 ((z) = 0.09), Tripp ct al. 2000 {{z) S3 0.24), Buries 
& Tytler 1996 ({z) = 0.9), and Simcoc, Sargent, & Rauch 2002 ((z) sa 2.5). 
The f2(Ei on ) values in column (5) are recalculated from the published data 
so that these values correspond to a flat cosmological model with Q m = 0.3, 
Qa = 0.7, and h = 0.65. Column (2) gives the fraction of the critical mass 
density contributed by the element; (4) the average redshift of the Lya 
systems from which column (5) is obtained; (5) the fraction of the critical 
mass density contributed by the ions of the element; and (6) the ionization 
fraction /(Ei on ) = f2(Ej on )/r2p(E) estimated from the P inventory. 



